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A cosmological model of dark energy interacting with dark matter and another 
general component of the universe is considered. The evolution equations for coinci- 
dence parameters r and s, which represent the ratios between the dark energy and 
the matter and the other cosmic fluid, respectively, are analyzed in terms of the sta- 
bility of stationary solutions. The obtained general results allow to shed some light 
on the coincidence problem and in the equations of state of the three interacting 
fluids, due to the constraints imposes by the stability of the solutions. We found 
that for an interaction proportional to the sum of the DE density and the third fluid 
density, the hypothetical fluid must have positive pressure, which leads naturally to 
a cosmological scenario with radiation, unparticle or even some form of warm DM 
as the third interacting fluid. 
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I. INTRODUCTION 

The existence of a dark component with an exotic equation of state, i. e., with a ratio 
w = p/p negative and close to —1, which drives an accelerated expansion is consistent with 
the luminosity distance as a function of redshift of distant supernovae [lj], the structure 
formation (LSS) {2] and the cosmic microwave background (CMB) j^] 

The cosmic observations show that densities of dark energy and dark matter are of the 
same order today To solve this coincidence problem [4J (or why we are accelerating in the 
current epoch due that the vacuum and dust energy density are of the same order today 
?) it is assumed an evolving dark energy field with a non-gravitational interaction with 
matter (5] (decay of dark energy to matter). 

Although the main topic of investigation have been centered in the interactions in the dark 
sector, it is physically reasonable and even expected from a theoretical point of view, that 
dark components can interact with other fluids of the universe. For example DE interacting 
with neutrinos was investigated in |fij], and decaying into the fermion fields, in A more 
general scenario was considered in [8|, in which DE is interacting with neutrinos and DM. 

Inspired in these previous investigation we have recently formulated a effective model 

where DE is decaying in DM and another hypothetical fluid 9], which we shall denote here 

and after by DX. In the framework of the holographic DE, and using the Hubble radius as 

infrared cutoff, we have shown that our scenario leads naturally, for a flat universe, to a 

more suitable approach to the cosmic coincidence problem in which the ratio between the 

energy densities of DM and DE, r, can be variable during the cosmic evolution. Our model 

has been discussed has a possible approach to solve the triple coincidence problem in 10]. 

In this work it was assumed that the third fluid is radiation. Although matter and radiation 

are almost non-interacting fluids, since the decoupling era, nevertheless they could interact 
I 

with DE. In was investigated the dynamical behavior when DE is coupling to DM and 
unparticle in the flat FRW cosmology. 

The aim of this paper is to investigate further the model of DE interacting with DM and 



another hypothetical fluid. Despite of the interesting results found in 10, III] when this third 
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fluid was specifically identified with radiation and unparticle, we continue assuming only 
that this fluid has an equation of state with uj constant. We expect, studying the stationary 
solutions of the evolution equation for r and s, the radio between the energy densities of 
DE and DX, to found suitable constraints for the equation of state of the unknown fluid. 
Since we are introducing a third fluid in the interactions in the dark sector, we generalize the 
coupling terms that have been already considered in the literature, although and interaction 



term which include the energy density of unparticle was discussed in [11] . In this paper we 
choose three different coupling terms in order to investigate the dynamical behavior of the 
models of DE interacting with DM and a third unknown fluid. 

Our paper is organized as follows. In section II we present the model for a universe filled 
with dark matter, dark energy and another fluid. We shall impose that the interacting term 
5 and 5' which appears in the conservation equations are different. In section III we study 
the stationary solutions and its stability for the equations of evolution of the ratios between 
the DM and DE and the third fluid and DE. We presents the constraints on the equation 
of state for DE and the hypothetical third fluid. In section IV we discuss our results. 



II. INTERACTING DARK ENERGY 



In the following we modeled the universe made of CDM with a density p m and a dark 
energy component, pr>. We will assume that the the dark matter component is interacting 
with the dark energy component, so their continuity equations take the form 

p D + 3H{l+u D )p D = -5', (1) 

p m + 3H{1 + uo m ) p m = 5. (2) 

where 5 and 8' are the interactions terms, which are different in order to include the scenario 
in which the mutual interaction between the two principal components of the universe leads 
to some loss in other forms of cosmic constituents. We assume that the interactions terms 
are of the form 5 = 3HT, where T = T(p D , p m , p x )- Since we will study an universe with 
dark energy decaying into dark matter, we have from the beginning 5 > and 5' > 0. If we 
denote this other component by px, its corresponding continuity equation is given by 

p x + 3H(l+u x )p x = 5' -5 = 3HU, (3) 
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where II = V — V. We are taking about in this case that dark energy decay into dark matter 
(or viceversa, depending on the sign of S) and other component. The sourced Friedmann 
equation is then given by 

3k 

3H 2 = p D + Pm + p x - — . (4) 

a 2 

Since we are interested in the three cosmic coincidence, we do not we assume as it was 
done in [9], that Pd ^ Px and p m 3> px- In order to study the evolution of the densities 
of these three fluids we construct the differential equations for the coincidence parameters 
t = Pm/ Pd and s = Px/pd- These equations takes the following expressions 



, r r p m p D 

r 



H H \p m Pd 
and 



(5) 



, S S ( p X PD \ ( . 

We assume that the interaction term takes the general form 

r = XiPd + A 2 p m + \ 3 px = (Ai + A 2 r + \ 3 s) p D . (7) 

Consistently, for the general form given in Eq.(j7]), we can parameterized II = V — F through- 
out suitable constants AJ, A2 , A3, in the following form 

n = Xlp D + \lp m + Xlpx = (XI + Kr + Xls) p D . (8) 

Notice that the form of interaction considered here is a generalization of the cases previously 
investigated, since the interactions taken account consider only functions of the dark sector 
nteractions that are linear combinations of the dark sector densities have been 



densities. 

studied in 1_2{], when the interaction is only between the dark components. 

Using Eqs.flU) and (TSJ) in the equations (jSJ) and (ED, with the interaction terms, Y and IT 
given by Eqs. ([7j) and jSJ), respectively, we obtain the evolution equations for the parameters 
r and s 



and 



r> = 3r [(Ax + A 2 r + A 3 s)(l + 1/r) + (XJ + X*r + X^s) +u D - u m } , (9) 



s' = 3s [Ai + A 2 r + X 3 s + (1 + l/s)(X n 1 + X^r + X n 3 s) + lu d - lu x ] ■ (10) 
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III. STATIONARY SOLUTIONS 

Before to look for the stationary solutions of Eqs. and (jSj) , let us briefly discuss the what 
we assume for the equations of state of the three cosmic fluids that are under interaction. At 
this stage, we only consider that one fluid, with an energy density pu and equation of state 
u)d, is decaying into the other two fluids. Although, it is reasonable to take u m = from 
the beginning, since we are thinking in the dark matter fluid, we shall postpone this election 
until to obtain the constraint derived from the study of the stationary solutions of Eqs.flS]) 
and (El) and its stability. For simplicity, the three equations of state are taken constant. 

Since we are looking for stationary solutions of Eqs.flHJ) and (jEJ) we set r' = s' = 0, 
obtaining a systems of algebraic equations in terms of the the variables r and s, with the 
parameters Ai, A 2 , A 3 , X±, A£ , A3, u D , u X - In order to study how suitable interactions leads to 
these three cosmic fluids to have stationary solutions we begin with the simplest cases. 



A. Case T = Xpo and II = X* pn 

In this case we are choosing A2 = A3 = A£ = A3 = 0. The interactions terms are 
proportional to the dark energy density. This type of interaction has been investigated 
in 11, Q, Q. The condition r' = s' = leads to the algebraic systems 

f(r, s)\ r=rs>s=Sa = A(l + r s ) + X n r s + (u D - u m )r s = 0. (11) 

g(r, s)\ r= r atS=Sa = Xs s + A 7r (l + s s ) + (u D - lo x )s s = 0. (12) 

Notice that the interaction assumed gives two simple linear equation non acopled in the 
variables r s and s s , which are the stationary solutions given by 

A 



A + A 71 " + u D - uj r , 

and 



(13) 



(14) 



A + X 71 + LJ D - UJ X 

Since r s and s s are positive quantities, the denominators of the above both equations must 
be negative, so be obtain the following inequalities 

uo D < uo m - (A + A"), (15) 
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and 



u x > u D + {\ + A 77 ). 



(16) 



In order to study the stability of the found solutions we shall evaluate the eigenvalues of 
the matrix M 



df(r,s) aj(r.s) 
M = I dr 8s 

dg(r.s) dg(r,s) 
dr ds 



(17) 



whose elements are evaluate at the critical point (r s , s s ). From the equation for the eigen- 
values, det[M. — 7/1], and since df(r, s)/ds = dg(r, s)/dr = 0, we obtain 



df(r, s) 
dr 

The eigenvalues are then 

m = 

and 

V2 = 



- V 



\r — r s ,s — s s 



df(r, s) 



dg(r, s) 
ds 



-v 



0. 



dr 



dg{r, s) 
ds 



A + A 7 " + LU D - UJ r , 



A + y + UJ D - ijJ X 



(18) 



(19) 



(20) 



The condition of stability, r\\ < and r/ 2 < is the same just contained in Eqs. f|T5|) and 
( JIB"]) . So, if for a given u>d, uj m , ux, A, A 71 there are positive stationary solutions, they 
are also stable. For the case of dark matter with negligible pressure, i.e., u m = 0, Eq. ffTS]) 
implies that the dark energy must necessarily have an equation of state with up < 0. A 
rough estimation of the value of the sum A + A 71 can be obtained from the equation for the 
acceleration 

- = ~(l + 3 W )p, (21) 
a o 

where p = pc(l + r + s) and to = {u)£> + oj m r + uj x s)/1 + r + s. If an accelerated phase is 
demanded we obtain the following inequality 



0Jd < -- [1 + (1 + 3cu m )r + (1 + 3u x )s} . 

In terms of the parameters Qd and fix (for u m = 0) we obtain 

1 



wo < 



3tt 



D 



(22) 



(23) 
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If it is assumed that fix -C 1, as we expect for today for any other fluid different from the 
dark sector. In order to obtain a rough estimation of the term A + X n , we equate the right 
hand side of the expressions ( fT5l) and ( 1231) . and taking Qr> = 0.7 we obtain 

A + A" ~ 0.5, (24) 

so for io n < —0.5 the three interacting fluids leads to stationary and stable solutions for r 
and s. From the inequalities (TT5l) and (jl~6l) we obtain that ujx ^ oj m . So for uj m = the 
third fluid could be a normal fluid or even an exotic fluid. Nevertheless, from the expression 
for the ratio ^, given by 

S = t( 1_ ^ + a + aO ' (25) 

we can conclude, assuming ^ « 1, that an scenario with u x > is more suitable taking 
A >> A 77 . From Eq. (fT3l) and taking r s m 0.3/0.7, up ~ —1, we obtains for the realistic 
case A >> A 7 " that A ~ 0.3. 



B. Case T = Xp m and II = X 7T p ri 



In this case we are choosing Ai = A3 = A^ = A3 = 0. The interactions terms are 
proportional only to the dark matter density. This type of interaction was investigated for 



models of interacting phantom dark energy with dark matter 
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15, 0,0 and also in Q, 



19| | . Observational constraints on A for this type of interaction have been investigated 
The condition r' — s' — leads to the algebraic systems 



f(r,s)\ 



r=r s ,s=s s 



r s 2 (A + A") + r s (\ + tu D - uo m ) = 0. 



(26) 



g(r, s) \ r=rsjS=Ss = s s r s (X + A 71 ") + X w r s + (u D - u x )s s = 0. (27) 

In this case the interaction assumed gives two coupled non linear equations in the variables 
r s and s s . The Eq. ( 1261) has the following solution different from zero , 

A + ujd — w m 
rs ~ A + A- • 

Imposing the condition r s > we obtain the constraint 



(28) 



X + LU D -uj m <0. (29) 
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Introducing the value for r s , given by Eq. (l28l) . in the Eq. (l27I) yields 

A + u m - u x V A + A- J K J 

Using the constrain given in Eq. (1291) in the expression for s s we obtain that s s > implies 

- A + uj m - uj x < 0. (31) 
From the equation for the eigenvalues, det\Wl — rjl], and since df(r, s)/ds = 0, we obtain 



df(r, s) 



dr 



- V 



r — r s ,s — s s 



dg(r, s) 
ds 



-7) 



\r — r s ,s — s s 



0. 



The eigenvalues are then 



and 



df(r, s) 
dr 



dg(r, s] 
ds 



\r= r s ,s — s 



2r s {\ + W T ){\ + uj D -u) TI 



r s (X + X n ) + lo d - uj x 



(32) 



(33) 



(34) 



Notice that the condition of stability, r/i < and r^ 2 < gives us, for r/i < the following 
constraint 



A + UJo — > 0, 



(35) 



which can not be allowed if Eq.( l29l) is satisfied. Then for this type of interaction it is 
not possible to have stable solutions for three cosmic interacting fluids. Nevertheless, this 
situation is a consequence of chose A > and A 71 " > from the beginning. It is straightforward 
to prove that for A + A 71 " < and A 71 " > 0, or equivalently A < and |A| > \X*\, the fixed 
point of the system are stable. The constraints for the equations of state are the following 



and 



A + uj d - oj m > 0, 



A + u D - ujx < 0, 



(36) 



(37) 



which for u m = yields A + uj d > and A + uj x > 0. Since A < the above conditions 
implies a third fluid with u>x > and also a DE with ujd > 0. So in the approach of DE 
interacting with two fluids, this kind of coupling give stable solutions but a cosmic evolution 
without acceleration. 
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C. Case r = Xp D (l + s) and 11 = X"p D (l + s) 

In this case we have taken Ai = A3, A2 = A£ = , and A^ = A3. The interaction is, in this 
case, proportional to po + Px- A coupling term which include a different fluid from those 



of the dark sector was already introduced in [U], but throughout expressions like pnPx and 



PmPxi where DX was identified with unparticle. The condition r' = s' = leads to the 
following algebraic system 

f(r, s)| r=rs , s=Ss = A(l + r s )(l + s s ) + AV S (1 + s a ) + {u D - u m )r s = 0. (38) 

g{r, s)\ r=rsjS=Sa = Xs s (l + s s ) + A 7r (l + s s ) 2 + (u D - u x )s s = 0. (39) 

Solving first the Eq. fl3"9l . which is a second grade equation of the form x 2 + Bx + C = 0, 
where the coefficients B and C are given by 

A + 2A- + u, -^ C= _Z_. (40) 
A + A 71 " A + A w 



The solutions of Eq. (1391) has the form 



Since s s is a positive and real number we need to impose the two constraints B < and 
B 2 > AC. The first one implies that A + 2A 7r + up — uox < and the second one, A + 2A 77 + 



oj d — uj x > —2a/ A 7r (A + A 71 ), which gives the following range for uj d — uj x 

- ((A + 2A") + 2a/>(A + A-)) < uo D - uo x < -(A + 2A 7r ). (42) 



Introducing the two solutions of Eq. (j39j) . which we denote by s s+ and in Eq. fl38|) . we 
obtain two solutions for r s , r s+ and r s _, given by 

A A 



(43) 



J- 1 Ss-\- J-i^s — 

Since r s > Eq. fj43l) gives the following constraint 

(A + A") (1 + s s ±) +u D -u m <0, (44) 

so u D must satisfy 

lj d <^ m -(A + A 7T )(l + s s± ). (45) 
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The conditions which are necessary to hold in order to have positive solutions for r s ± and 
s s ± are then the inequalities given by Eqs. fj42l and f|44l) 

As in the case A, where the interaction is proportional only to the DE density, if uj m = 
the above inequality implies that the dark energy must necessarily have an equation of state 
with ujd < 0. A rough estimation for the range of the values that the parameters A and 
A 71 " can take, may also be done for this case using Eq. (1231) with fix ^ 1 and Qd = 0.7. 
Equating the right hand side of the expressions fj4"5l) and fl23|) and since 1 + s s ± > 1 



A + X* < 0.5. 



(46) 



Let us to study now the conditions imposed on the found solutions if we require stability. 
Evaluating the elements of the matrix M at the critical points (r s+ ,s s+ ) and (r s _,s s _), we 
obtain from the equation for the eigenvalues, det[lSA — rjl], and since dg(r, s)/dr = 0, that 



df(r, s) 
Or 



\r=r„±,a=s B ± 



dg(r, s) 
ds 



0. 



The eigenvalues are then 

'df(r, s) 



Or 



r— r ,s- 



and 



V'2 



dg{r, s) 
ds 



(X + X 7T )(l + S s± )+UJ D -LJ ri 



(A + 2A 7r )(l + s s ±) + Xs s± + uo D - uo x . 



(47) 



(48) 



(49) 



|r = r s± ,s = s s± 

The condition of stability, rji < and 772 < 0, implies that the following constraints must 
hold 



and 



uo D - uo m < -(A + A^)(l + s s ±) 



uo D -uo x <- [(A + 2A 7r )(l + s s ±) + Xs s± ] 



(50) 



(51) 



Notice that the constraint given by Eq. fl50l) is the same obtained in Eq. fj4*4|) . Nevertheless, 
we need to look for the range of ujd — which can accommodate the constraint given by 
Eq.( l42l and Eq.( l5TI) . Choosing Eq.(l511 as the constraint for the upper limit of ujd — u>x, 
the upper limit indicated in Eq. (l4"2"|) is also satisfied. We can impose the condition 



((A + 2X") + 2v/A-(A + A-)) < - [(A + 2A 7r )(l + s s± ) + Xs s± ] 



(52) 
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which leads to the following condition for s a ± 




A + \* 



(53) 



Since it is physically reasonable assume s s ± <C 1 for a late time evolution, 

It is straightforward to check from the expressions Eq. (j501) and Eq. flSTj) that, indepen- 
dently of the critical point considered, the equation of state of the non decaying fluids satisfy 



If uo m = 0, which is the equation of state for the dark matter fluid, the above result indicates 
us that the exigency of stability for the stationary solutions of the evolution equations (Q 
and (jTOj) . imposes an unknown interacting fluid with non null pressure. 



In the present investigation we have considered a cosmological scenario where the dark 
energy is decaying into the dark matter and another component of the universe, which we 
do not identify explicitly. We have assumed that this three fluids have an equation of state 
with u constant. We have choose three different coupling terms, analyzing the stationary 
solutions of the evolution equation for parameters r and s. 

When the coupling is proportional to the dark energy only, we have found that the 
conditions for the stationary solutions be positive are the same of those to be stable. For 
dark matter with negligible pressure we obtain that the dark energy must necessarily have 
an equation of state with ojd < 0. 

When the coupling is only proportional to the DM energy density, since in this case 
it is not possible to obtain stable solutions for the three interacting fluids, if A > and 
A 71 " > 0, which guaranty that DE decaying in the other fluids. Nevertheless, relaxing this 
condition and taking A + A 71 " < and A 71 " > we obtain stationary solutions which are stable. 
Notice that this physically correspond to a DM decaying in DE and in the third fluid. It is 
interesting to mention that the usual case of DE interacting only with DM, which have been 



discussed for this coupling in [20||, showed that the data slightly favored a DM decaying in 
DE and ujd < —1. Unfortunately in our approach, we obtain that the unknown fluid and DE 
have positive pressure, leading to a decelerated expansion. As in the case of DE interacting 
only with DM, there is no stable solutions for DE with negative pressure. 



(54) 



IV. DISCUSSION 
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In the third case studied, the coupling consider a different fluid from those of the dark 
sector, taking a term proportional to the sum of the DE density and the third fluid density. 
We have found two fixed points and from the constraints derived from the condition of 
stability we have obtain that ujx > w m , which means an interacting third fluid with positive 
pressure. This type of coupling can then accommodate an scenario with radiation 10], 
unparticle 11] or even some form of warm DM 2l| as the third interacting fluid. 
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